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Abstract 

We investigate the ground-state magnetic order of the spin-1/2 J1-J2 Heisenberg model 
on the square lattice with ferromagnetic nearest-neighbor exchange J± < and frustrat- 
ing antiferromagnetic next-nearest neighbor exchange J2 > 0. We use the coupled-cluster 
method to high orders of approximation and Lanczos exact diagonalization of finite lattices 
of up to N = 40 sites in order to calculate the ground-state energy, the spin-spin correlation 
functions, and the magnetic order parameter. We find that the transition point at which the 
ferromagnetic ground state disappears is given by J^ 1 = 0.393|Ji| (exact diagonalization) 
and J^ 1 = 0.394|Ji| (coupled-cluster method). We compare our results for ferromagnetic J\ 
with established results for the spin-1/2 J1-J2 Heisenberg model with antiferromagnetic J\. 
We find that both models (i.e., ferro- and antiferromagnetic J\) behave similarly for large 
J2, although significant differences between them are observed for J2/I-/1I < 0.6. Although 
the semiclassical collinear magnetic long-range order breaks down at J| 2 ~ O.6J1 for anti- 
ferromagnetic Ji, we do not find a similar breakdown of this kind of long-range order until 
J2 ~ 0.4|Ji| for the model with ferromagnetic J\. Unlike the case for antiferromagnetic J\, 
if an intermediate disordered phase does occur between the phases exhibiting semiclassical 
collinear stripe order and ferromagnetic order for ferromagnetic J\ then it is likely to be over 
a very small range below J2 ~ 0.4| J\\. 

PACS codes: 
75.10.Jm Quantized spin models 

75.45. +j Macroscopic quantum phenomena in magnetic systems 

75.10.Kt Quantum spin liquids, valence bond phases and related phenomena 

1 Introduction 

The spin-1/2 Heisenberg antiferromagnet with nearest-neighbor (NN) J\ and frustrating next- 
nearest-neighbor (NNN) J2 coupling (J1-J2 model) on the square lattice has attracted much at- 
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tention during the last twenty years (see, e.g., Refs. [TH2T)] and references therein). This model may 
serve as a canonical model to study the interplay of frustration effects and quantum fluctuations 
as well as quantum phase transitions driven by frustration. 
The corresponding Hamiltonian reads 

H = Ji^Sf Sj + J 2 ^2 Si ■ sj, (1) 

where the sums over and run over all NN and NNN pairs, respectively, counting each 

bond once. 

Recently, several quasi-two-dimensional magnetic materials with a ferromagnetic (FM) NN 
coupling Ji < and an antiferromagnetic (AFM) NNN coupling J 2 > have been studied 
experimentally, e.g., Pb 2 VO(P0 4 ) 2 (2TH2I], (CuCl)LaNb 2 7 [25], SrZnVO(P0 4 ) 2 [MIEHIET] , and 
BaCdVO(P04) 2 |23l[26j|28] . Due to a quite large AFM coupling J 2 these materials are driven out 
of the FM phase. These experimental studies have stimulated a series of theoretical investigations 
of the ground state (GS) and thermodynamic properties of the FM Ji- J 2 model, i.e. the model 
with FM NN exchange J x and frustrating AFM NNN exchange J 2 [2DH5T] . It was found that the 
FM GS for the spin-1/2 model breaks down at J 2 = J 2 cl « 0.4| J x | [2HHIIHI31EIIEE]. Note that 
for the classical model (s — > oo) the corresponding transition point is at J 2 = 0.5| J\\. Moreover, 
for sufficiently large J 2 > J% 2 a magnetically long-range ordered collinear stripe GS appears. 
According to Refs. [30 | l3T | l35 l l36] this second critical frustration strength is J% 2 ~ 0.6 1 Ji|, and both 
magnetically ordered phases are separated by a magnetically disordered phase with enhanced 
nematic correlations. Note that this collinear stripe GS phase, as well as the values of J% 1 and J^ 2 
reported in the literature, are similar to the observation for the corresponding AFM model, i.e. 
the model with AFM J x [Dia[6l[71|9l[nHl3l[l5H20] . 

As we know from the investigation of the AFM J\- J 2 model the study of the GS phases in the 
strong frustration regime, i.e. J 2 ~ 0.5Ji, is one of the hardest problems in the field of frustrated 
quantum magnetism. In the meantime, scores of papers have dealt with this problem. On the 
other hand, for the FM Ji-J 2 model so far only a few studies exist. Hence, further investigations 
using alternative approaches are highly desirable to confirm or to query the existing results. 

Due to frustration, highly efficient quantum Monte-Carlo codes, such as the stochastic series 
expansion suffer from the minus sign problem. Therefore, many other approximate methods, e.g. 
the Green's function method [IQ1C321EI], the series expansion [H1II21IIZ], the Schwinger boson 
approach [7], variational techniques [9], the functional renormalization group technique [2D] as 
well as the projected entangled pair states method [18] were used to study the GS phases of the 
model with AFM J x . 

In the present paper we calculate the GS energy, spin-spin correlation functions, and the 
magnetic order parameter of the collinear stripe order using two completely different methods, 
namely a large-scale exact diagonalization (ED) of finite lattices up to N = 40 sites (see Sec. [2]) 
and the coupled-cluster method (CCM) (see Sec. [3]). Both methods have been successfully applied 
to the AFM Ji-J 2 model, see e.g. Refs. [M1CE9] for the ED method and Refs. [HUME] for the 
CCM. 

Henceforth we set | J\\ = 1 if not stated otherwise explicitly. 

2 Results of Lanczos exact diagonalization for finite square 
lattices 

The paper of Schulz and co-workers [1], presenting results for finite lattices up to iV = 36 sites for 
the first time has set a benchmark for ED studies of quantum Heisenberg magnets. Frequently, such 
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Table 1: Singlet GS energy per site E /N and square of order parameter M^(0,7r), cf. Eq. fl2]), 
for the Ji- J 2 model with J\ = — 1 on finite square lattices of N = 32, 36 and 40 sites. 





iV = 32 


N = 32 


N = 36 


N = 36 


N = 40 


iV = 40 


J2 


E /N 


M 2 N (0,ir) 


E /N 


M 2 n (0,tt) 


E /N 


M&(0,7T) 


0.40 


-0.31284253 


0.11167341 


-0.31431081 


0.08925821 


-0.31127637 


n 


10687524 


0.45 


-0.34234578 


0.12363200 


-0.34259033 


0.11439947 


-0.34062128 





11332095 


0.50 


-0.37406800 


0.12617129 


-0.37369533 


0.11854042 


-0.37174835 





11581048 


0.55 


-0.40673497 


0.12730259 


-0.40577614 


0.12049011 


-0.40374930 





11719238 


0.60 


-0.43994975 


0.12793318 


-0.43841728 


0.12170653 


-0.43628038 





11810286 


0.65 


-0.47352394 


0.12833506 


-0.47143169 


0.12256200 


-0.46917022 





11876286 


0.70 


-0.50735237 


0.12861459 


-0.50471467 


0.12320526 


-0.50231993 





11927048 


0.80 


-0.57553353 


0.12898259 


-0.57184472 


0.12411474 


-0.56916790 





12000667 


0.90 


-0.64418634 


0.12921694 


-0.63949153 


0.12472863 


-0.63652042 





12051714 


1.00 


-0.71315309 


0.12937976 


-0.70748822 


0.12516916 


-0.70421577 





12089143 



numerical exact results are also used to test new approximate methods, see, e.g. Refs. [Tt3 | l3"9~ | l4"0"] . 
Due to the progress in the computer hardware and the increased efficiency in programming, very 
recently the GS and low-lying excitations of the unfrustrated (i.e., J2 = 0) [H] as well as the 
frustrated spin-1/2 HAFM [19] and of a spin-1/2 Heisenberg model with ring exchange [12] have 
been calculated by the Lanczos algorithm for a square lattice with N = 40 sites. The largest 
two-dimensional quantum spin model for which the GS has been calculated so far is the spin-1/2 
HAFM on the star lattice with N = 42 sites [13]. Note, however, that for sectors of S z > which 
are relevant for finite magnetic fields much larger system sizes can be treated by ED, see, e.g., 



In order to analyze GS magnetic ordering, the spin-spin correlation functions and an appro- 
priate order parameter are important quantities. Following Schulz et a!., [1] we use here the 
Q-dependent susceptibility (square of order parameter) defined as 



The relevant order parameter for the collinear stripe magnetic LRO present at large J2 is Mjy(Q) 
at the magnetic wave vectors Qi = (vr,0) or Q 2 = (0,7r). In Table [1] we give the singlet GS 
energies per site as well as the order parameters Mjy(0, ir) for finite lattices of N = 32, 36, and 40 
sites. In a second Table [2] we present the spin-spin correlation functions (so • Sr) for the largest 
lattice considered. There are in total 11 different correlation functions for N = 40 that are given 
in Table [2] for the same data points as in Table [TJ Periodic boundary conditions are imposed for 
all of the lattices employed here. A graphical representation of the finite lattices can be found in 
Refs. [1] and [19]. The exact data, which is provided here in some detail in the tables, might be 
used as benchmark data for approximate methods applied to this modelJl] 

The corresponding data for the energies, the square of order parameters and for some selected 
correlation functions are also displayed in more detail in Figs. [U [21 and [3j By way of comparison, 
we also present corresponding curves for the AFM model (Ji = +1) [19J. We recall that the GS 

lr The spin-spin correlation functions for N = 20, 32, and 36 can be provided upon request. 



Ref. gl]. 
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Table 2: Spin-spin correlation functions (s • Sr) for the J1-J2 with J\ = — 1 model on the finite 
square lattice of N = 40 sites. 





(s ■ s R ) 


(s ■ s R ) 


(s ' s R ) 


(s • s R ) 


(s • s R ) 


(s ■ Sr) 


J2 


R=(1,0) 


R= (2,0) 


R= (3,0) 


R=(3,l) 


R = (2,l) 


R=(l,l) 


0.40 


0.044832 


0.183933 


0.012219 


-0.194442 


-0.040218 


-0.277017 


0.45 


0.033144 


0.196362 


0.014787 


-0.199197 


-0.033852 


-0.304812 


0.50 


0.027648 


0.201318 


0.014973 


-0.200208 


-0.029595 


-0.316452 


0.55 


0.024204 


0.204078 


0.014523 


-0.200667 


-0.026460 


-0.323037 


0.60 


0.021753 


0.205887 


0.013884 


-0.201012 


-0.024012 


-0.327312 


0.65 


0.019875 


0.207189 


0.013203 


-0.201327 


-0.022023 


-0.330324 


0.70 


0.018363 


0.208185 


0.012537 


-0.201621 


-0.020367 


-0.332568 


0.80 


0.016044 


0.209616 


0.011328 


-0.202149 


-0.017745 


-0.335676 


0.90 


0.014313 


0.210600 


0.010293 


-0.202593 


-0.015750 


-0.337719 


1.00 


0.012954 


0.211320 


0.009417 


-0.202962 


-0.014172 


-0.339156 





(s • Sr) 


(s • Sr) 


(s ' Sr) 


(s ■ Sr) 


(s • Sr) 


J2 


R=(l,2) 


R=(l,3) 


R= (2,2) 


R= (2,3) 


R = (4, -2) 


0.40 


-0.027909 


-0.177408 


0.160527 


-0.007815 


0.154329 


0.45 


-0.023862 


-0.183867 


0.176142 


-0.003522 


0.166329 


0.50 


-0.021057 


-0.186297 


0.181923 


-0.002268 


0.169644 


0.55 


-0.018936 


-0.187719 


0.185043 


-0.001671 


0.171225 


0.60 


-0.017247 


-0.188748 


0.187071 


-0.001323 


0.172227 


0.65 


-0.015858 


-0.189567 


0.188526 


-0.001098 


0.172968 


0.70 


-0.014691 


-0.190248 


0.189642 


-0.000936 


0.173562 


0.80 


-0.012825 


-0.191325 


0.191262 


-0.000723 


0.174483 


0.90 


-0.011394 


-0.192135 


0.192387 


-0.000588 


0.175173 


1.00 


-0.010260 


-0.192765 


0.193215 


-0.000498 


0.175707 
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Figure 1: GS energy per site E (N)/N calculated by ED for finite lattices of N = 32,36,40 sites 
as well as the extrapolated value eo (N — > oo), cf. Eq. flH]). For comparison we show also E (N)/N 
for the AFM model (Ji = +1) with N = 32, 36, 40 [T9]. The blue line shows the energy of the FM 
eigenstate. 

is identical for the FM and the AFM model on the classical level for J 2 > 0.5. For the spin- 1/2 
model, due to quantum fluctuations, the two models behave differently. Only in the limit J 2 — > oo 
the sign of J\ becomes irrelevant for the quantum GS. Our numerical data for the GS energy, 
the order parameter, and the spin-spin correlation functions demonstrate that for J 2 > 0.8 the 
different models behave qualitatively very similarly. Note, however, that due the different sign in 
the NN exchange J\, also the NN spin-spin correlation (s ■ Sr.), R = (1, 0), has a different sign for 
the models. In the region 0.8 > J 2 > 0.6 the differences between both models become more and 
more pronounced. For J 2 < 0.6 both models behave completely different. From many previous 
studies [THT9] we know that J 2 « 0.6 is the point where for the AFM model the semiclassical 
collinear stripe LRO gives way for a magnetically disordered phase. It is obvious that at J 2 « 0.6 
the order parameter Mjy(0, 7r) of the AFM model changes drastically (cf. Fig. [21 whereas M^(0, tt) 
remains almost constant up to about J 2 « 0.45 for the FM model, i.e. there is no indication for a 
breakdown of semiclassical magnetic LRO around J 2 ~ 0.6 for J\ = — 1. Only for J 2 < 0.45 there 
is a noticeable decrease in the order parameter M^(0,7r) until the transition point J 2 l between 
the singlet GS and the fully polarized FM GS. Note, however, that M^(0,7r) remains finite until 
J 2 l . The different behavior of both models can also be seen in the GS energy and the spin-spin 
correlation functions, cf. Figs. [1] and [3j From Fig. [1] it is evident that the finite size-effects in 
the GS energy Eq of the FM model are small. In fact, the difference in Eq between N = 36 and 
N = 40 is less than 1% in the whole region shown in Fig. [TJ Hence, a finite-size extrapolation 
should give accurate results for the GS energy (see below). 



5 





0.14 




1 ? 

U. 1 c 




0.1 




0.08 


cT 






0.06 








0.02 









32, Jh=-1 
36, J 1= -1 
40, J 1= -1 
40, J 1= +1 



0.5 0.6 



J 2 /J 1 



Figure 2: Square of order parameters M^(0,7r), see Eq. (j2J), for iV = 32,36, and 40 (Ji = 
For comparison we also present a corresponding curve for the AFM model (J\ = +1) for iV 



0.3 



A 

DC 
CO 

o 

CO 
V 



S— B— B— E] 




Figure 3: Selected spin-spin correlation functions (s • Sr) calculated with Lanczos ED for N 
(thick lines: J\ = —1, thin lines J\ = +1). 
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Figure 4: Extrapolated order parameters m s , see Eq. (Hj), for the FM model (Ji = —1) and the 
AFM model (J x = +1). 

Following the lines of Refs. jHEEH] we use now the data for N = 20,32,36,40 to perform a 
finite-size extrapolation. Note that we do not include the data for the lattice with N = 16 sites, 
since this lattice exhibits an anomalous behavior jl]. The finite-size extrapolation rules for the 
considered two-dimensional model are well known [151 HO]- The finite-size behavior of the GS 
energy is given by 

E (N) _ const. 

— - e o - -jpn + ■ ■ ■ • (3) 

The extrapolated energy eo is shown in Fig. [TJ We use the intersection point between the ex- 
trapolated singlet GS energy e and the energy of the fully polarized FM GS to determine the 
transition point between both GS phases to J^ 1 = 0.393. 

The order parameter for the collinear stripe LRO in the thermodynamic limit is defined as [41TT9] 
m s = y/8 liniAr-s.oo Mjv(0, 7r). The finite-size behavior of Mn(0, 7r) is given by 

, ^9 , \ 1 9 const. 

il4(vr,0) = -m2 + -^ + ... • (4) 

However, it is evident from Fig. [2] that the finite-size behavior of M^(tt,0) becomes irregular at 
J 2 = 0.44, since we have Mf (7r, 0) > Mf 6 (7r,0) for Jf < J 2 < 0.44. Hence, the extrapolation 
of M^(7r,0) becomes inconclusive in that region. The results for m s for J 2 > 0.44 are shown in 
Fig. HI For comparison we show also the extrapolated order parameter for the AFM model |19j . 
Although m s becomes zero at J 2 2 ~ O.66J1 for the AFM model, we do not find an indication for 
a magnetically disordered phase in the region above J 2 = 0.44 for the FM model. However, we 
cannot exclude such a GS phase in the very small parameter region located between J 2 = 0.393 
and J 2 = 0.44 from our ED data. 
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Figure 5: The GS energy per spin as function of J2 obtained by CCM-LSUBn with n = 4, 6, 8, 10 
and extrapolated values in the limit n — > 00 using the extrapolation scheme e(n) = ao + ai(l/n) 2 + 
a2(l/ri) 4 . For the sake of comparison we show also the extrapolated GS energy for the AFM 
model [16]. The blue line shows the energy of the FM eigenstate. 
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3 Results of the coupled-cluster method (CCM) for the 
infinite square lattice 

For the sake of brevity, we will illustrate only some important aspects of the CCM. The interested 
reader can find more details concerning the application of the CCM on the AFM J1-J2 in Refs. [HI 
[T3HT6]. For more general information on the methodology of the CCM, see, e.g., Refs. (I7H39]. 

First we notice that the CCM approach automatically yields results in the thermodynamic 
limit N — > 00. The starting point for a CCM calculation is the choice of a normalized reference 
(or model) state |<&). We then define a set of mutually commuting multispin creation operators 
Cj with respect to this state, which are themselves defined over a complete set of many-body 
configurations /. For the system under consideration here we choose one of the two degenerate 
classical collinear stripe states as the reference state. Next, we perform a rotation of the local 
axis of the spins such that all spins in the reference state align along the negative z axis. In 
the rotated coordinate frame the reference state reads |$) = | \) \ \) \ \) . . . , and we can treat 
each site equivalently. The corresponding multispin creation operators then can be written as 
Cj = sj, sjsj, sjsjsj, ■ • • , where the indices i, j, k, . . . denote arbitrary lattice sites. 

The CCM parametrizations of the ket- and bra- GS's are given by 

H\V) = E\V); (y\H = E(&\; 
|M/) = e 5 |$), S = Y,SiCf; 

(V\ = (<5>\Se- s , S=l + J2^i C T- ( 5 ) 

We wish to determine the correlation coefficients Si and Si for the correlation operators S and 
S. By using the Schrodinger equation, H\ty) = E\ty), we can write the GS energy as E = 
(<&\e~ s He s \<&). The magnetic order parameter is given by 

1 - 

i=l 

where si is expressed in the rotated coordinate system and s = 1/2 is the spin quantum number. 
To find the ket-state and bra-state correlation coefficients we have to solve the so-called CCM ket- 
and bra-state equations given by 

(<S>\Cye- s He s \<S>) =0, V/ ^ 0, (7) 
($\Se- s [H, C/]e 5 |$) = 0, V/ ^ 0. (8) 

Each ket- or bra-state equation belongs to a certain creation operator Cj' = sj, sjsj, sjsjsj, ■ ■ • , 
i.e. it corresponds to a certain set (configuration) of lattice sites i, j, k, . . . . 

For the considered quantum many-body model we have to use approximations in order to 
truncate the expansion of S and S. We use the well established LSUBn scheme P^HTBlHTHoT] in 
which all multispin correlations in the correlation operators S and S over all distinct locales on 
the lattice defined by n or fewer contiguous sites are taken into account. For instance, within the 
LSUB4 approximation, multispin creation operators of one, two, three or four spins distributed on 
arbitrary clusters of four contiguous lattice sites are included. The number of these fundamental 
configurations can be reduced by using lattice symmetry and conservation laws. In the highest 
order of approximation considered here, namely the LSUB10 approximation, we have finally 45825 
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Figure 6: Selected spin-spin correlation functions (so ■ sr) calculated with CCM within LSUB6 
approximation (thick lines: J\ = — 1, thin lines J± = +1). 



fundamental configurations for the collinear stripe reference state, yielding 45825 coupled nonlinear 
equations which have to be solved numerically. 

The LSUBn approximation becomes exact for n — > oo, and so we can improve our results 
by extrapolating the "raw" LSUBn data to n — > oo. There is ample experience regarding how 
one should extrapolate the GS energy eo(n) and the magnetic order parameter m s (n). For the 
GS energy per spin e (n) = a + a\(l / n) 2 + a 2 (l / n) 4 is a reasonable well-tested extrapolation 
ansatz p^1 - [T6ll48T[5T] . An appropriate extrapolation rule for the magnetic order parameter is 
m s (n) = b + 6i(l/n) 1/2 + b 2 (l/n) 3/2 pUJIfi] . Moreover, we know from Refs. [T4Tfi~6] that the 
lowest level of approximation called the LSUB2 approximation conforms poorly to these rules. 
Hence, as in previous calculations P3HT6], we exclude LSUB2 data from the extrapolations. 

Note that starting from large J 2 the solution of the CCM-LSUB10 equations can be traced 
until J 2 = 0.402, only. For J 2 < 0.402 no CCM-LSUB10 solutions with respect to the collinear 
stripe reference state could be found. This could be related to the complexity of the set of 45825 
coupled equations and the fact that the fully polarized FM state is a low-lying excitation very 
close to the true GS for J 2 ~ 0.4. However, we cannot exclude the possibility that this problem 
might also be an indication that the collinear stripe order is not a good description for the true 
GS for J 2 values very close to the transition to the FM GS. 

In Fig. [5] we show the GS energy per site for the LSUB6 and LSUB10 approximations, as well 
as the extrapolated GS energy (n — > oo). The numerical CCM-LSUBn data demonstrate, that the 
difference between the various LSUBn data is very small (see Fig. [5]). As a result, the extrapolated 
GS energy almost coincides with LSUB10 data (in fact, the difference is less than 0.4%). Moreover, 
the difference between the extrapolated CCM energy and the finite-size extrapolated ED energy 
is less than 1% in the whole parameter region. Hence, we conclude that the GS energy data 
presented in this paper for the FM J\-J 2 model are very accurate. We again use the intersection 
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Figure 7: Magnetic order parameter m s versus J 2 obtained by CCM-LSUBn with n = 4, 6, 8, 10 
and its extrapolated values (n — > 00) using the extrapolation scheme m s {n) = b + ^(l/n) 1 / 2 + 
b 2 (1/n) 3 / 2 . By way of comparison we show also the extrapolated collinear stripe order parameter 
(n -> 00) for the AFM model (J x = +1). 

point of the energies in order to find an estimate for the transition point Jf 1 to the FM GS. We 
find that J 2 l = 0.395 for the LSUB6 approximation and this is shown in Fig. [5j We find also 
that J 2 l = 0.394 and J 2 l = 0.398 for the LSUB8 and LSUB4 approximations, respectively (not 
shown in Fig. [5]). As the difference between the LUSB6 and LSUB8 results for J 2 l is already very 
small we may choose the LSUB8 value, J^ 1 = 0.394, as the CCM estimate for the transition point. 
This CCM critical value is in excellent agreement with the ED estimate J 2 l = 0.393 (see above). 
We also compare the GS energies for the FM and the AFM models in Fig. [5j As for the ED 
method, the deviation in energies between the two models (i.e., negative and positive J\) becomes 
substantial at around J 2 = 0.6. 

We also present the CCM results for spin-spin correlation functions for the FM as well as the 
AFM model in Fig. [6j where we have shown the same correlation functions as in Fig. [3j Recall 
that one of the two degenerate collinear stripe states has to be chosen as CCM reference state. As 
a result the NN correlation functions (s • sr), R = (1, 0), and (s • sr), R = (0, 1), are different, 
whereas, e.g., all NNN correlation functions are the same. Hence, we show in Fig. [6] the averaged 
NN correlation function. For the FM model the agreement of the CCM correlation functions with 
the ED data is excellent, see Figs. [3] and [61 For the AFM model there is a noticable difference 
for J 2 < 0.59. In particular, the steep decrease in the correlation functions at J 2 = 0.59 present 
in the ED data for the AFM model, cf. Ref. [19], is not observed in the CCM LSUB6 data. 
For J 2 < 0.59 most likely no magnetic LRO exists for the AFM model [16J. Therefore it is not 
surprising that the CCM solution in a finite order of LSUBn approximation based on the collinear 
stripe reference state does not povide quantitatively precise results for the correlation functions 
inside the magnetically disordered phase. 

Last but not least, we present results for the collinear stripe order parameter m s calculated 
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by the CCM via Eq. ([H]) in Fig. [7J For the order parameter, the CCM-LSUBn results depend 
on the approximation level n more strongly than for the GS energy. Hence, it is more important 
to extrapolate them in the limit n — > oo. By comparing the correlation functions and the order 
parameter of the FM model with the AFM model, it is again obvious that both models behave 
similarly for larger J 2 . However, there are significant differences for J 2 < 0.6 between both models. 
Although the CCM result for the AFM model clearly shows the breakdown of the semiclassical 
collinear magnetic LRO around J 2 = 0.6 (the CCM estimate for the critical frustration is Jf 2 ~ 
0.59 [IS]), the spin-spin correlation functions and the order parameter for the FM model stay 
almost constant over the whole parameter region shown in Fig. [7J Again, we obtain good agreement 
of the order parameter between the ED data and the CCM data in a wide parameter range, cf. 
Figs. [Hand [7J Finally, we come to a similar conclusion as in Sec. |2]that we do not find an indication 
for the breakdown of the collinear stripe magnetic LRO until J 2 ~ 0.4 for the FM model. 

4 Summary 

In this paper we have presented results on the GS magnetic ordering of the J\- J 2 spin-1/2 Heisen- 
berg magnet with FM NN exchange J\ = —1 on the square lattice obtained by Lanczos ED 
as well as the CCM. In agreement with previous studies [2lMJIll3"4"] l3~7 ] l3"8"] we obtain values for 
the transition point J| x at which the fully polarized FM GS (present for small J 2 ) breaks down, 
J 2 cl = 0.393 1 Ji | (ED) and J 2 cl = 0.394| J x | (CCM). 

In contrast to previous studies [2"9"H3"T| |3"4"] we do not find indications for the breakdown of 
the semiclassical collinear magnetic LRO (present for large J 2 ) at about J 2 = 0.6 1 J\\, rather this 
magnetic LRO seems to be stable till J 2 ~ 0.4|Ji|. 

Although our results are in favour of semiclassical magnetic GS long-range order for a wide 
range of frustrating exchange J 2 , the low-lying excitations and, as a consequence, the low-temperature 
thermodynamics might be strongly influenced by frustration. This becomes increasingly important 
when approaching the transition to the FM GS at J 2 X ~ 0.4 1 Ji|, see e.g. Ref. [37]. For J 2 > Jf 1 
the FM multiplet becomes a low-lying excitation, and so an additional low-energy scale will appear 
leading, e.g., to an extra low-temperature peak in the specific heat, see, e.g., Refs. [29 | l52 | [53] . 

Stimulated by the recent experimental activities and the related search for new square-lattice 
materials near the quantum critical point of the Ji-J 2 model, the issue of low-lying excitations 
and their relevance for the thermodynamic properties should be investigated in future studies in 
more detail. 
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